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Abstract 

Using the two-particle self-consistent approach and cluster perturbation theory for the two-dimensional t-t'-t"-U 
Hubbard model, we discuss weak- and strong-coupling mechanisms for the pseudogap observed in recent angle 
resolved photoemission spectroscopy on electron-doped cuprates. In the case of the strong-coupling mechanism, 
which is more relevant near half-filling, the pseudogap can be mainly driven by short range correlations near the 
Mott insulator. In the vicinity of optimal doping, where weak-coupling physics is more relevant, large antiferro- 
magnetic correlation lengths, seen in neutron measurements, are the origin of the pseudogap. The t — J model is 
not applicable in the latter case. 
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Angle resolved photoemission spectroscopy [1] 
(ARPES) provides deep insight into the nature 
of high-temperature superconductors. In partic- 
ular, it has revealed the failure of Fermi liquid 
theory to describe single-particle excitations in 
these systems. Contrary to the quasiparticle con- 
cept of Fermi liquid theory, certain segments of 
the would-be Fermi surface are almost gapped. 
This is the so-called pseudogap phenomenon. In 
particular, recent ARPES measurements [2] on 
Nd2-a;Ce x Cu04 have shown that, in contrast to 
the hole-doped cuprates, lightly electron-doped 
(e-d) ones have a large spectral weight near (n, 0). 
With further doping towards optimal doping, 
spectral weight also appears around the zone di- 



agonals, leaving hot spots (regions with large scat- 
tering) in between (ir, 0) and (tt/2, tt/2) where the 
non-interacting Fermi surface intersects the anti- 
ferromagnetic (AFM) zone boundary. Theoretical 
explanation of these experimental data is still an 
open question. 

In the present paper we discuss two approaches 
for the pseudogap in e-d cuprates: a strong and 
a weak-coupling one. In the case of the strong- 
coupling mechanism [3], which is more relevant 
near half-filling, the pseudogap is mainly driven 
by short-range correlations near the Mott insu- 
lator. In the vicinity of optimal doping, where a 
weak-coupling mechanism is more appropriate, 
large AFM correlation lengths are the origin of 
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the pseudogap [4]. 

We use two different methods, the two-particle 
self-consistent (TPSC) approach [5,6] and cluster 
perturbation theory [7] (CPT), for the single-band 
Hubbard model on a square lattice with a repul- 
sive local interaction U and nearest t, next-nearest 
t' and third-nearest t" neighbour hoppings. The 
former method is based on a self-consistent de- 
termination of the irreducible vertices that enter 
dynamical susceptibilities (spin-spin and density- 
density). This is done by enforcing the Pauli prin- 
ciple, conservation laws for spin and charge fluctu- 
ations, and important sum rules. These results are 
then used to obtain an improved approximation for 
the single-particle self-energy. TPSC has been ex- 
tensively checked against Quantum Monte Carlo 
simulations [5,6,8]. The CPT approach is based on 
exact diagonalization of finite clusters that are cou- 
pled through strong-coupling perturbation theory. 
The CPT results were calculated on 4 x 4 clus- 
ters. Since TPSC is valid in the weak to interme- 
diate coupling regime, while CPT is more reliable 
for intermediate to strong coupling, we can study 
all values of U to gain insight into two mechanisms 
for the pseudogap in e-d cuprates. 

Weak-coupling pseudogap: In the vicinity 
of optimal doping, the absence of zero energy 
excitations around (tt/2,it/2) is a general result 
of strong-coupling calculations [3,9] for t-t'-t"-U 
and t-t'-t"-J models. This is consistent even with 
weak-coupling calculations [10,11] for the Hub- 
bard model when the coupling is increased towards 
the value of the bandwidth. On the other hand, 
several papers [3,4,10,12,13] have shown that the 
spectral weight near (tt/2,it/2), as seen in exper- 
iment [2], shows up at optimal doping when the 
Hubbard coupling U is somewhat smaller than the 
bandwidth. Therefore, one concludes that strong- 
coupling physics is not relevant for the pseudogap 
in c-d cuprates near optimal doping, and that a 
weak to intermediate coupling mechanism is ap- 
propriate. In other words, the t-J model does not 
describe the physics of e-d cuprates near optimal 
doping. Further evidence for this is the disap- 
pearance at x = 0.15 of the lower Hubbard band 
that is present at x — 0.04 in ARPES data [2] on 
Nd2- K Ce a; Cu04. One can derive the t-J model 
from the Hubbard model only when both upper 



and lower Hubbard bands are well defined. 

At weak coupling, large AFM correlation 
lengths, seen in neutron measurements [14], are 
the driving force for the pseudogap [4] . The physi- 
cal mechanism is that electrons on a planar lattice 
suffer scattering by AFM fluctuations which have 
large phase space in two dimensions [5,15]. Thus, 
those quasiparticles in regions of the Fermi sur- 
face that can be connected by the AFM vector 
(so called hot spots) do not survive scattering by 
these strong AFM fluctuations. This theory [4] 
explains in detail the ARPES results mentioned 
above, as well as the temperature dependent corre- 
lation length measured by neutron scattering [14]. 
Ref. [4] also makes a few predictions for ongoing 
experiments: (a) The ARPES pseudogap found 
at low temperatures should be seen even in the 
paramagnetic phase up to the temperature when 
the AFM correlation length becomes smaller than 
the single-particle thermal de Broglie wavelength 
£th — frvF /irksT '. The corresponding pseudogap 
temperature T* is close to that found in optical 
experiments [16]. (b) For T < T* and for T some- 
what larger than T*, the characteristic spin fluc- 
tuation energy in neutron scattering experiments 
is smaller than the thermal energy (renormalized 
classical regime) and the spin fluctuations are 
overdamped near T* . 

In the spirit of Ref. [4] we also compare the 
renormalized Fermi velocities at optimal doping 
n = 1.15 along the (it, 0)-(7T, it) direction and along 
the zone diagonal with the corresponding ARPES 
data [17] on Ndi.gsCeo.isCuO^ Using the experi- 
mental renormalization factors and bare Fermi ve- 
locities [17], the experimental renormalized Fermi 
velocities are 3.31 x 10 5 m/s and 3.09 x 10 5 m/s 
along the zone diagonal and along the (it, 0)-(tt, n) 
direction, respectively. The corresponding renor- 
malized Fermi velocities obtained by TPSC are 
3.27 x 10 5 m/s and 2.49 x 10 5 m/s, respectively. 
The agreement is very good, particularly along the 
diagonal direction. The bare Fermi velocities are 
renormalized in TPSC by roughly a factor of two. 

In Fig. 1 we present the single-particle spectral 
weight obtained by TPSC for filling n — 1.175 and 
temperature T = i/50 for the same values of U 
and band parameters as in Ref. [4]. Here, A(k,u>) 
is multiplied by the Fermi-Dirac distribution func- 
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Fig. 1. Energy distribution curves A<(k.,u>) = A(k, ui)f(u>) 
along the Fermi surface shown in the inset for n = 1.175, 
U = 5.75t. The lines are shifted by a constant for clar- 
ity. Band parameters are t 1 = —0.1754, t" = 0.05t while 
T = t/50. 



t'=-0.175t, t"=0.05t 



t'=-0.3t, t"=0.2t 




C75t, n=1.175 
(d) 



CVI 








00 




6 










CPT 


° 


■ 


0% 



tion f(oj). One can see that A < (kp,oj) is peaked 
at zero energy near (tt, 0) and {tt/2, tt/2), and is 
shifted away from the Fermi energy (pseudogaped) 
towards higher binding energies at hot spots where 
the Fermi surface intersects the AFM Brillouin 
zone boundary. The calculated Fermi surface plot 
.A<(k, 0) that corresponds to this case (Fig. 2(d)) 
shows, at hot spots, the zero energy suppression 
of spectral weight by AFM fluctuations. Thus, we 
predict that the pseudogap induced by AFM fluc- 
tuations should be experimentally seen up to 18% 
electron doping, and should disappear for larger 
dopings [4]. Note however that 18% doping in our 
calculations may be equivalent to 15% doping in 
reduced samples if we take the point of view of 
Ref. [14] that the role of reduction can be mod- 
eled as a Ax w 0.03 shift with respect to nomi- 
nal Ce concentration. Then our data for the Fermi 
surface plot and spectral function at n = 1.175 
should be compared with the ARPES ones for re- 
duced samples near x — 0.15 Ce doping. The semi- 
quantitative agreement of Ref. [4] is preserved since 
Fermi surface plots and energy distribution curves 
for n = 1.15 (see Figs. 1,2 of Ref. [4]) and n = 
1.175 (shown here) look similar. However, in better 
agreement with experiment, the theoretical pseu- 
dogap feature for 15% doping is slightly more pro- 
nounced and extends over a broader region in k 
space than that for 17.5% doping. For comparison 
of our correlation lengths with neutron measure- 



Fig. 2. (color) Fermi surface plots A<(k, 0) in the first 
quadrant of the Brillouin zone by two methods: TPSC and 
CPT, T = (0,0), M = (tt.tt). The TPSC plot in part (a) 
and (b) are at T = 0.05t, that in part (d) is at T = 0.02i 
while the CPT plot in part (c) is at T = 0. 

ments data on 15% doped as grown samples we 
should use n = 1.15 as was done in Ref. [4]. Note 
however that Fig. 3 of Ref. [16] may lead us to argue 
that doping corresponds to nominal Ce concentra- 
tion only in reduced samples. The role of reduction 
and Cc alloying on doping remains to be clarified. 

As shown in the Fermi surface plot obtained by 
TPSC [Fig. 2(a)], for t' = -0.175t, t" = 0.054, 
strong AFM fluctuations at 12.5% doping cause 
the suppression of spectral weight not only at hot 
spots, but also along a large segment of the Fermi 
surface near (tt/2, 7r/2). This is consistent with 
strong-coupling calculations [3] . The AFM correla- 
tion length is about 40 lattice spacings for this plot, 
which is larger than ^ t h, and the spin susceptibil- 
ity at (tt, tt) is much larger than the noninteracting 
one. By contrast, for large values of |4'| = 0.34 and 
t" = 0.24, one can see in Fig. 2(b) the Fermi liq- 
uid like (uniform) distribution of spectral weight 
in k space that results from the strong suppres- 
sion of AFM fluctuations by frustration. Here the 
AFM correlation length is about two lattice spac- 
ings, and the spin susceptibility at an incommensu- 
rate wave vector (jr — 5, tt) is just a few times larger 
than the noninteracting one. This again shows that 
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Fig. 3. Single-particle spectral weight obtained by CPT at 
n = 1.125 as a function of energy lu in units of t for wave 
vectors along the high symmetry directions shown in the 
inset. Band parameters are t' = — 0.3t, t" = 0.2t. 

large AFM correlation lengths are the driving force 
for the pseudogap in the weak-coupling case. 

Strong-coupling pseudogap: An additional 
type of physics becomes relevant in the strong- 
coupling regime, which we argued is more appro- 
priate when electron doping is decreased towards 
half-filling. In this limit, short range correlations 
are sufficient to create a pseudogap even in the 
absence of long AFM correlation lengths [3]. As 
shown in the Fermi surface plot obtained by CPT 
[Fig. 2(c)], for large values of \t'\ = 0.3i and t" = 
0.2£ (enough to frustrate antiferromagnetism in 
TPSC) the pseudogap does occur at 12.5% dop- 
ing around the nodal direction as a result of strong 
coupling, whereas it is absent at weak coupling for 
the same values of if and t" [see Fig. 2(b)]. We can- 
not rule out the existence of long AFM correlation 
lengths at strong coupling with CPT because the 
finite cluster size precludes such correlations. 

The physical mechanism for the pseudogap here 
is that those quasiparticles in regions of the Fermi 
surface that are connected to other such regions 
by wave vectors that have a broad spread of ra- 
dius 8 around (tt, tt) suffer strong scattering by 
short range correlations [3]. Thus, the pseudogap 
driven by the Mott physics with short range cor- 
relations occurs around zero energy and only in 
these regions of the Fermi surface (hot spots). In 
contrast, the Mott gap occurs for all wave vectors 
and is not tied to zero energy. This is illustrated on 
Fig. 3, which shows energy dispersion curves ob- 



tained by CPT at n = 1.125 for wave vectors along 
the high symmetry directions shown in the inset. 
The range of frequencies away from zero energy 
where A(k, u>) — for all wave vectors is the Mott 
gap. At finite electron doping it always opens up 
at negative energies when the Hubbard coupling U 
is sufficiently large, and the chemical potential lies 
in the upper Hubbard band. The range of frequen- 
cies around zero energy where spectral weight is 
suppressed only for some wave vectors along T-M 
and X-H directions is the pseudogap. This suppres- 
sion of spectral weight at zero energy is also seen 
in the Fermi surface plot for the same parameters 
[Fig. 2(c)]. Note that the longer range AFM cor- 
relations, which are observed in experiment [14], 
would probably only reinforce the strong-coupling 
mechanism that already exists in the presence of 
short-range correlations. One may speculate that 
at strong coupling the condition £ > £ t h becomes 
£ > a (with a the lattice spacing), which is easier 
to satisfy. This would connect the weak and strong 
coupling regimes in a continuous manner. 
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